Dipole spectrum structure of non-resonant non-pertubative driven two-level atoms 
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We analize the dipole spectrum of a two-level atom excited by a non-resonant intense monochro- 
matic field, under the electric dipole approximation and beyond the rotating wave approximation. 
We show that the apparently complex spectral structure can be completely described by two families: 
harmonic frequencies of the driving field and field-induced nonlinear fluorescence. Our formulation 
of the problem provides quantitative laws for the most relevant spectral features: harmonic ratios 
and phases, non-perturbative Stark shift, and frequency limits of the harmonic plateau. In partic- 
ular, we demonstrate the locking of the harmonic phases at the wings of the plateau opening the 
possibility of ultra-short pulse generation through harmonic filtering. 
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I. INTRODUCTION 

Two-level atoms are essential blocks for the under- 
standing of basic processes in quantum physics. In par- 
ticular they are the keystone of quantum optics as the 
harmonic oscillator is the basic model for classical op- 
tics. Their simplicity brings the possibility for analyti- 
cal developments which convey physical interpretations 
not attainable with more exact approaches. Moreover, 
two-level models have good quantitative accuracy in the 
description of laser-matter interactions near resonance at 
moderate laser intensities, where the rotating wave ap- 
proximation is applied pQ. For the case of strong fields, 
where a large number of atomic or molecular transitions 
are simultaneously involved and ionization is present, the 
two-level model does not provide a complete description, 
although it can still be used to study the role of bound 
state transitions in high-order harmonic generation [2J. 
Nevertheless, under particular conditions some of these 
complex systems are dominated by two-level transitions 
even for the strong driving case (for instance, the molec- 
ular Hydrogen ion at large internuclear distances [3J). 
Two-level transitions have also been associated to har- 
monic generation in molecules involving a rescattering 
center different from the electron's parent ion [3] and con- 
nected with charge-resonant states in odd-charge molec- 
ular ions [5]. Nowadays the recent trends in OP- CPA 
(Optical Parametric Chirped-Pulse Amplification) tech- 
niques to produce high-power mid-infrared laser radia- 
tion [5] renew the interest of few-level systems in the 
off-resonant strong-coupling regime, since in this spectral 
region few atomic transitions can be reasonably isolated 
from the rest [3 [5] • In addition, the forthcoming devel- 
opment of extreme intensity laser sources may permit to 
directly address the excitation of dipolar two-level tran- 
sitions in nuclei [5]. 

Two-level atoms interacting with intense laser light 
show a significant non-linear behavior dominating the 



complex temporal dynamics. This dynamics is translated 
to the atom dipole evolution and, therefore, to the struc- 
ture of the emitted field. The radiation spectra under 
such circumstances shows a plateau structure, in which 
harmonics with similar intensities extend to high frequen- 
cies |10j . This is a universal behavior, present in other 
models of light-matter interaction (also in classical anhar- 
monic oscillators) and the signature of the failure of the 
perturbative regime. The physics beneath this structure 
is well understood for the case of atoms and molecules in 
strong electromagnetic ionizing fields [TT]. In this case, 
the relevant process responsible for the plateau consists in 
the ionization of the electrons, their acceleration through 
the field-induced quiver motion, followed by the recolli- 
sion with the parent ion with the release of the electron's 
kinetic energy into electromagnetic radiation. On the 
other hand, for a two-level system, the plateau emerges 
from transitions between adiabatic states as was pointed 
out in (T2J [T3] (see also the dicussion in the last section 
of this paper). 

The study of coherently driven two-level systems be- 
yond the perturbative limit has a fundamental reference 
in |14j . with special focus on the near resonant condi- 
tions. In Ref. [15] . a quantitative description of the har- 
monic spectra of strongly driven two-level systems was 
reported by means of a Floquet-based theory. This ap- 
proach resulted in a closed iterative analytical formula- 
tion, based on a continued-fraction method, that pro- 
vided for exact ratios between the harmonic intensities. 
Later on, this approach was used to study the harmonic 
enhancement connected with field-induced multiphoton 
resonances 16J. Other continued fraction approaches 
have also been used to derive approximated analytical 
expressions for the plateau extension [T7]- The non- 
harmonic peaks also present in the radiative spectrum 
of the electric dipole have been studied for the case of 
large frequencies and / or very high amplitudes of the driv- 
ing fields [IH1 EI]- In these studies, the non- harmonic 
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resonances are identified as satellite peaks around each 
harmonic, separated by the transition frequency between 
the quasienergy states. Two groups of peaks were iden- 
tified (hyperRaman lines of Stokes and antiStokes type) 
according to their relative position with respect to the 
harmonic peaks. Also, these authors derived analytical 
formulas for the quasi-energies. In the present paper we 
use the exact approach developed in [TS] to address the 
complete description of the dipole spectrum in the gen- 
eral case. We will show that the apparent complexity of 
the spectra can be reduced to two contributions, namely, 
harmonic radiation and laser-induced fluorescence. This 
latter family accounts for the emergence of the mentioned 
satellite structure (Stokes and antiStokes lines) in the 
case of strong driving and/or in the limit of high driving 
frequencies. In particular, our approach gives a simple 
explanation of the asymmetry between the intensities of 
the Stokes and antiStokes satellite peaks. On the other 
hand, we obtain a polynomial equation whose solution 
gives the position of the Stark-shifted transition energy 
to, in principle, arbitrary accuracy. Restricting to the 
lowest degree of this equation, we derive an analytical 
formula for this later quantity. We test it against the 
exact numerical integration of the two-level Schrodinger 
equation, and show the convergence to the result of [18] 
(in the limit of validity considered in this reference). We 
also derive approximated formulas for the extension of 
the harmonic plateau. Our formula for the cut-off en- 
ergy agrees with |17j and converges to |19j in the limit 
of Rabi frequencies much greater than the transition fre- 
quency. On the other hand, our formula for the plateau's 
on-set frequency converges to |17j in the low frequency 
limit and gives better comparison with the exact results 
for higher driving frequencies. Finally, from the phase 
of the higher-order harmonics, we show their locking and 
demonstrate numerically the possibility of producing at- 
tosecond pulses. 

II. THEORY 

We begin with the standard formulation for a two-level 
system driven by a monochromatic field, out of the rotat- 
ing wave approximation. The dynamics of the probability 
amplitudes are given by: 

ij f ai(t) = -l[xe- iut +X*e iut ]a 2 (t), (1) 

ij t a 2 {t) = cu a 2 (t) - \ [ X e-^ + X *e i " t ] o x (t) , (2) 

where labels 1 and 2 stand for the lower and upper level 
respectively, ujq is the atomic transition frequency, u> is 
the monochromatic field frequency, and x — £o/V^ i s 
the Rabi frequency, where £q is the electromagnetic am- 
plitude and fi the dipole matrix element. The electro- 
magnetic amplitude is defined by a linearly polarized 
monochromatic field of the form E(i) = E cos(o;< + </>), 
being Eq and 4> the amplitude and the phase of the field 



respectively, with £ — |E |e -1 ^, and the dipole matrix 
element as [i = (l\ez\2), being a real quantity and as- 
suming z the field polarization direction. We define the 
Bloch variables u and v as the real and imaginary parts 
of ai{t)a2{t), respectively, and the population inversion 
as w(t) = \a 2 (t)\ 2 -\a 1 (t)\ 2 . Then, from Eqs. (l) and {2} 
it can be derived: 

ju{t) = -uov(t) , (3) 

j t v{t) = uou(t) + \ [xe- 1 ^ + xV" 4 ] w(t) , (4) 

jw{t) = -2[xe-^+x*e^]v(t). (5) 

The atomic dipole is defined as d(t) = 2fj,u(t). Figure [I] 
shows the spectral content |w(w)| 2 computed numerically 
from Eqs. (|3|-([5|. To be specific, in this paper we have 
chosen the driving field phase as <fi = — vr/2. To reduce 
the effect of the abrupt envelope, we have considered a 
sinus squared field turn-on of two cycles, followed by 16 
cycles of constant amplitude. The dipole spectrum has 
been calculated through the Fourier transform of u(t) 
in these later 16 cycles. For small field amplitudes, fig- 
ure [l|a), the spectrum is dominated by two peaks, one 
located at the driving field frequency uj (Rayleigh scat- 
tering) and another at the transition frequency (field- 
induced fluorescence). As the field amplitude increases, 
figure Qb), the complexity of the dipole spectrum in- 
creases with the appearance of harmonic frequencies of 
the driving field surrounded by satellite peaks, see in- 
set of figure [ljb) , whose position varies with the driving 
field intensity. These satellite structures are referred as 
hyperRaman lines in [IS] [05] . We should note that the 
relative intensity of the two satellite peaks around each 
harmonic varies with the driving field parameters, the 
satellite structure being asymmetric in the general case. 
Also, note that the harmonic structure shows the emer- 
gence of a plateau region with peaks of similar intensities. 
Mid-Infrared radiation can provide a suitable scenario in 
order to match the parameters considered in Fig. [l]^b). 
For example, the 2s-3p transition in the Hydrogen atom is 
close to 5 photon resonance using a 3 micron wavelength 
laser. In this case, the current technology is almost ready, 
see for example Ref. _J5]. 

Combining Eqs. (J3j) and Q we obtain 

^u(t) = -f [ X e~ iut +x*e iut ] w(t) , (6) 

which describes the atomic dipole evolution. As dis- 
cussed in [T5], the dipole dynamics corresponds to the 
harmonic oscilator in the small coupling limit (w(t) ~ 
— 1) and, therefore, eq. ^ provides a connection be- 
tween the fundamental atomic models of classical and 
quantum optics. On the other hand, combining Eqs. ^ 
and ([5]), we have 

^0 = A U . e -... + A v»']|„ W . (7) 
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FIG. 1: Dipole spectrum, |u(tj)| 2 , of the driven two-level 
atom, (a) spectrum at low field amplitude (Ixl/k-'o = 0.005 
and oj/cjo = 0.1), (b) spectrum at higher field amplitude 
(|x|/^o = 1-5 and uj/ujq = 0.2). Inset: detail of the dipole 
spectrum in (b) showing the structure of satellite peaks 
around the harmonics. 



an intermediate expression that we will use in the follow- 
ing. We now consider the spectral decompositions of u 
and w, 



u(t) = I u(a)e~ lat da , 

J — OO 

w{t) = / w(a)e- lat da. 



(8) 
(9) 



Note that we are now considering Fourier components 
of arbitrary frequency, instead of the harmonic series in 



[To] , restricted to integer multiples of the laser frequency. 
The spectral forms associated with Eqs. (|6| and ^ are 



u(a) 



w(a) 



wo 



cjq — a 2 w{a — uj) uJq — a 2 w(a + uj) 



UJ 



a — uj + a + uj 

X u(a — uj) + x u(a + uj) 

a a 



(10) 
(11) 



respectively. 



III. THE EXISTENCE OF TWO SPECTRAL 
FAMILIES 



The combination of the later expressions ( 10 1 and ( 11 ) 



leads to a set of (infinite) coupled algebraic equations 

$-{a)u(a-2uj)+e(a)u(a)+<I> + (a)u{a+2uj) = , (12) 
with 



9(a) = 2a 2 



(a 2 -uj 2 )(a 2 -uj 2 ) 



$±(a) 



(a ± 2uj)(a =F w) 



(13) 
(14) 



Note from these equations that the dipole spectrum is 
partitioned in a set of independent families of the sort 
{ao, ao ± 2w, ao i 4a;, • • • } (For instance, the set of har- 
monic peaks corresponds to a single family with ao = UJ )- 
This is a general result for two level systems interacting 
with a monochromatic field monochromatic fields, since 
up to this point we have not made any other approxima- 
tion. 

Each spectral family is described by the following tri- 
diagonal matrix, generated by Eq. (12), 



M(a) 



Q(a — 4u) $ + (a-4j) 

(a - 2uj) 9(q - 2uj) $+ (a - 2uj) 

$_(a) 6(a) $+(a) 

$_(a + 2a;) 9(a + 2a;) $+(q + 2cj) 
$_(a + 4a;) 9(a + 4w) 



(15) 



Since Eq. (12 1 is an homogeneous system of equations, 
the necessary condition for the existence of a family of 
spectral peaks is det(M) = 0, yielding a polynomial equa- 
tion (in principle of infinite degree), which gives exactly 
all the non-zero components of the dipole spectrum. The 



complete spectrum is, consequently, composed by all the 
families which are solutions of this equation. To find 
them we should address the general problem of finding 
all a so that det(M) = is fulfilled. 

Since any physical spectrum cannot have Fourier com- 
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ponents of arbitrarily large frequencies, all the spec- 
tral families must have some central region where the 
peaks have relevant intensities. In the following, we 
shall consider ao as a reference frequency belonging to 
this central region. Away of this region, the intensities 
of the spectral peaks should decay gradually as their 
frequency approaches ±00. This permits the trunca- 
tion of M at arbitrary large frequencies with almost 
perfect accuracy. Therefore, in the following, we shall 
consider M as truncated to dimension n x n (n odd) 
around the central frequency ao- In this case, each 
family of spectral peaks is composed by n components: 
{«0i £*o i 2cj, ao ±4w, • • • , ao ± (n — l)/2w}, each being a 
solution of det(M) = 0. Using definitions ( [l3| and (14), 
the truncated det [M (a)] = leads to an algebraic equa- 
tion of degree 4n in a. Therefore, the general spectrum 
will be composed at most by four families of peaks. Since 
the dipole is a real quantity, any solution for a implies 
a solution for —a (actually the associated Fourier com- 
ponent for —a is the complex conjugate of the one for 
a). Thus there are only two independent non- vanishing 
families composing the spectrum. As we have pointed 
out before, the harmonic field is already one of such fam- 
ilies, therefore there is only room for one more. We shall 
see later that this second family is associated with the 
natural frequency of the atom and, therefore, it is gener- 
ated by the non-linear mixing of the two-level transition 
frequency with the laser frequency. Note this is also a 
general result for monochromatic fields. 




FIG. 2: Dipole spectrum, |?t(cj)| 2 , plotted for negative and 
positive frequencies in two different cases (a) uj/ljo = 0.1, 
|xlM> = 0.5 , and (b) cj/lj = 0.2, |x|/w = 1.5. This later 
corresponds to the same case as in Fig. [ljb), including the 
inset. The two families composing the spectrum have been 
labeled with solid black triangles(harmonic family) and solid 
green boxes (fluorescence family). The complex conjugate of 
this later is labeled with solid red circles. 

Figure [2] shows the exact solution of the dipole spec- 
trum computed from Eqs. ([lj and ([2]). Despite it is 
always symmetric, for convenience we have plotted the 



negative as well as the positive frequency part of the 
spectrum. In plot (a), one can clearly identify the two 
families of peaks composing the spectrum: the harmonic 
family (filled black triangles) and the fluorescence family 
(filled green squares), the filled red circles correspond to 
the complex conjugate of this later. The harmonic family 
contains its own conjugate. The apparent complexity of 
the spectra at high field amplitudes, for instance the case 
shown in Fig. |T|(b) , raises when the fluorescent family ex- 
tends to negative frequencies, and its conjugate to posi- 
tive. In this case each harmonic appears as surrounded 
by two satellite peaks, one belonging to the fluorescence 
family and the other to its complex conjugate. Figure 
[2jb) shows the same case as[ljb), with the peaks labelled 
according to their correspondent family. As mentioned 
above, in general the intensity of the satellites around a 
particular harmonic is not the same, as the relative posi- 
tion of each satellite peak in the fluorescence family and 
its conjugate do not coincide. 



IV. STARK SHIFT 

The physical interpretation of ao in the fluorescence 
family can be found in the limit of low intensities % — > 0. 
As discussed before, in this limit the dipole evolution cor- 
responds to the classical harmonic oscillator. This later 
problem is described only by two spectral contributions, 
one peak at the laser frequency ui and other at the natu- 
ral frequency u>o of the oscillator. As the field increases, 
the two-level system reveals its non-linearity, and natu- 
ral and field frequency are mixed. Clearly, the harmonic 
spectrum rises from the mixture of the field frequency 
with itself, while the fluorescence family comes out of 
the mixture of the natural and field frequencies. There- 
fore, the central frequency ao of the fluorescence family 
is to be interpreted as the effective transition frequency 
of the two-level atom, which is AC Stark shifted from ujq 
as the field increases. 

Next we shall derive an approximated expression for 
the Stark shifted transition frequency. To do this, we 
evaluate the effective transition frequency ao by solving 
det(M) = 0. For an arbitrary large truncation of M, 
a limiting exact solution for this quantity can be found 
numerically. If we are interested, however, in closed ana- 
lytical estimations we should consider the regime of mod- 
erate coupling, for which the fluorescence family may be 
composed of several peaks but with a clear maximum 
at the central frequency ao- This situation allows us 
to neglect the influence of the peaks at the sides of a 
and to consider only the central term of matrix M (i.e. 
an 1 x 1 truncation). The effective transition given by 
det(M) = reduces to 6(a ) = 0, whose solutions are 



2|xl 



2|x| 2 ) 2 -4w\ 



(16) 
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FIG. 3: (a) Plot of formula (16 1 for uj/uiq = 0.1 versus the 
Rabi frequency, resulting from the truncation of M to 1 X 1. 
ao,± refer to the frequencies of the central peak in each of the 
two families composing the dipole spectrum. Since ao,- cor- 
responds to the harmonic family, it is practically insensitive 
to the intensity field . (b) Test of the analytical solution ao,+ 
against the exact numerical results from Eqs. |l]) and pj, 
for w/wo = 0.1. The numerical data shows a good agreement 
with the analytical formula for the Stark shift although for 
large intensities small deviations appear. Parts (c) and (d) of 
this figure show the results for ao,± resulting from the trunca- 
tion of M to 3 x 3. The results are obtained numerically from 
the equation det(M3 X 3) = and are compared with: (c) the 
exact numerical result also shown in part (b) of this figure, 
and (d) with the results for ao,- from the lxl truncation 
shown in part (a) of this figure. 



The stark shifted transition frequency corresponds to 
a o,+ (which has the correct limit ujq when |x| tends to 0), 
while ao,- converges to the laser frequency in the limit 
|x| tending to 0. Therefore, the two solutions in (16) 
represent the two possible spectral families, harmonics 



and fluorescence, which are restricted to single peaks as 
a result of the lxl truncation of M. In the limit of small 



field intensities, and for uj > 2lj, Eq. (16) is reduced to 
a Stark shift of \x\ 2 <jJo/(oJq — w 2 ), in correspondence with 
the form presented in [15] . 

Note that, even with the drastic truncation of M, the 
harmonic ao,_ is practically insensitive to the intensity 



field. In Fig. |3ja) we plot Eq. 1 16 ) for u/uq = 0.1 (taking 
only the positive frequencies), increasing the field from 
Ixl/^o = 0.01 to |xl/ w o = 0.5. The figure clearly shows 
that ao,- remains practically constant as the intensity 
field increases, at variance with oto,+> which clearly de- 
pends on the field intensity. The comparison of the ana- 
lytical solution ao.+ given by Eq. ( 16 ) with the numerical 
results of the fluorescence peak computing Eqs. ([I]) and 
^ give rise to a good match, see Fig. |3jb). As expected, 
the approximated formula obtained by ( |16[ ) is practically 
exact at low field intensities but gradually deviates, at 
large intensities, from the exact numerically calculated 
solution. 

As discussed before, a better accuracy can be found in- 
cluding the influence of the fluorescence peaks sorround- 
ing the central one at a. For instance, the nearest neigh- 
bors are included if we truncate M to a 3 x 3 matrix. In 
this case, det(M) = yields a polynomial of order 12 in 
a, that cannot be solved analytically. However, the nu- 
merical solution of the problem can be attained easily in 
this case, and also for higher truncation orders. In Fig. 
[3|c) we plot the Stark shift of the fluorescence transition 
and the points obtained by det(M 3x3 ) = 0. As expected, 
the accuracy of the results increases in this latter case 
compared with the order truncation n = 1. We depict 
in Fig. ao for the harmonic family in the n = 1 

and n = 3 truncation order. Note that the dependence 
of this quantity with the field strength is smaller in the 
latter the frequency of the harmonic should be 

independent of the intensity in the exact case (see the 
following section). 



V. THE HARMONIC FAMILY 

In order to analyze the harmonic spectrum, we con- 
sider the matrix (15) evaluated in ao = uj. First, we 
demonstrate that in this case the condition det(M) = 
is independent of the resonance frequency (w ) or the 



laser intensity (|x|)- Note from (14) that <3> + (a = oS) = 
and = — uS) = 0. These two zeros in M provide 

a supplementary symmetry that allows us to reduce the 
complexity of the matrix to the following structure 



M(a) 







$+(-3w) 
O(-Lu) $+(-w) 

$_(3u;) 



(17) 



where A (±w) are ,J -^ x submatrices (n is the trun- laser intensity. Now, starting from the new form (17) of 
cation of the matrix M), which depends on the laser fre- the truncated matrix M, we can calculate its determinant 
quency u, but also on the resonance frequency and the 



G 



as 



det[M] = det[A(oj)] ■ det[A(-u)] x 

[Q(u)Q(-cj) - = , 



(18) 



since [Q(o;)6(— ui) — $+(— cj)$_(w)] = 0, see equations 



(131 and (14). Therefore the location of the harmonic 
familiy is independent of the laser intensity. Let us now 
study the different aspects of the harmonic spectrum. 



A. Harmonic ratios and relative phases 



Coming back to Eq. 



tween neighboring harmonics as Z{a 
we have the following relation 



(12), and defining the ratio be- 
u(a) /u(a 



Z{a) 



$_(a) 



6(a) +$+(a)Z(a + 2uj) 



2w), 



(19) 



This recursive expansion is exact and, therefore, re- 
produces accurately the relative weights between the 
peaks of the harmonic and fluorescence spectral families 
[15] . While the ratios between the harmonics inside the 
plateau region form a complex sequence, it is not so in 
the regions outside the plateau, where the harmonics de- 
crease monotonically with the frequency, and the relative 
ratio Z(a) is a small quantity. In such regions we may 
approximate 



Z(a) 



$_(a) 



(20) 



The spectral region with frequencies below the plateau 
is characterized by the conditions a > 2a; and a < loq, 
yielding 0(a) > and $_ oc exp(— 2i<f)) . Consequently 
Z(a) oc — exp(— 2icj>) and the relative phase between con- 
secutive harmonics before the plateau is exp [—2i<fi + it]. 
On the other hand, the spectral region above the plateau 
is characterised by a arbitrarily large, therefore, 0(a) < 
and $_ oc exp(— 2i(f>). Consequently, the relative 
phase between consecutive harmonics after the plateau 
is exp(— 2i<j)). Inside the plateau region, the phase dis- 
tribution is in general more complex and can only be 
determined using the recurrent relation ( 19 ) analogously 
as it was done in [15] for the harmonic intensities. Fig- 
ure[4]shows the harmonic spectrum and phases computed 
numerically from Eqs. ([3| to ([5| in the strong coupling 
case (|x|/w = 4 and w/w = 0.2). The plateau region is 
shaded in grey (for a detailed discussion of the plateau's 
limits see below). As stated before, the field phase has 
been taken as cf) — — tt/2 and, therefore, the above dis- 
cussion predicts a relative phase between consecutive har- 
monics of 2tt before the plateau and tt after the plateau, 
in accordance with the numerical results shown in figure 
|1J Note that the above phase relations hold also approx- 
imately in the extreme parts of the plateau. This is a 
relevant aspect, as it implies that the harmonics near the 
plateau's cut-off are approximately phase locked. This 



spontaneous locking is also found in harmonic generation 
with ionizing systems, and implies the possibility of syn- 
thesis of attosecond pulses from the Fourier synthesis of 
the harmonics near the cut-off. This possibility is ana- 
lyzed in figure[5]in which we present the time evolution of 
the dipole resulting from the inverse Fourier transform of 
the spectral components with frequencies above 34u>, for 
the case shown in figure |4j The black filled curves show 
the squared envelope of the dipole, resulting in a train 
of pulses with durations well below the optical period. 
The irregularity (no periodicity) of the squared envelope 
of the dipole in Fig. [5] is because of the resonant family 
contribution. 




1 5 9 13 17 21 25 29 33 37 41 
Harmonic number 



FIG. 4: Dipole spectral intensities (a) and phases (b) for the 
strongly driven case with Ixl/^o — 4 and ui/ujq — 0.2. The 
harmonic peak family is highlighted with open circles. The 
shadowed box encloses the harmonic plateau structure. The 
limits of this box have been defined using the expressions ( 22 1 
and ( 24 1 for the plateau on-set and cut-off frequencies 



The 

arrow points to the plateau's on-set frequency in the weak 
driving limit, equal to ujq. 



B. Plateau's on-set and cut-off frequencies 

We can have an approximated idea for the plateau's 
extension by finding its limiting frequencies. To do this 
we will simply consider as a reasonable estimation the 
fact that the neighboring harmonics in the plateau have 
similar intensities: \u(a — 2uj)\ ~ \u(a)\ ~ \u(a + 2uj)\. 
The plateau on-set can be obtained using this condition 
together with the relative phase for the lower frequency 
harmonics derived above: exp(— 2i(f) + tt). For this case 
Eq. ([121 yields to 



$_(a)e 



+2i4 



0(a) + $ + (a)e 



-2i<j> 



, (21) 



7 



from Eq. (22 1 (corresponding to the lower limit of the 
shadowed box) is about 2w above ujq (pointed out by an 
arrow) . 
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FIG. 5: Black filled curve: detail of the inverse Fourier trans- 
form of the higher frequency part of the dipole spectrum 
shown in figure [4] Frequencies below 34cj are filtered out. 
Only the field envelope is represented. Red line: Sketch of 
the amplitude of the driving field at the same time interval. 



FIG. 6: Frequencies for the plateau's on-set (red line) and 
cut-off (green line) derived from Eqs. ( |22[ ) and ( |24| for the 
case w/cjo = 0.1. The visual estimation from the numerically 
computed spectra are plotted with black and blue points, we 
have added an error bar of ±2u as an error estimation of the 
method. 



which has a solution 
o 1 



a 



w, 



o 



+ lu 2 + J (u 2 - uj 2 ) ^ + 16w 2 |x| I 



(22) 



In the limit u> — > 0, the plateau's on-set is approximately 
w which is the estimation in |17) . 

The value of the frequency at the plateau's cut-off can 
be inferred imposing the phase condition for the higher 
frequency harmonics: exp(— 2i<fr). We have then 

$_(a)e+ 2 ^ + 9(a) + $+(a)e~ 2# = , (23) 

which has a solution 



w 2 



■4|x| : 



(24) 



in coincidence with [17] and also converges to [19] in the 
limit of Wo ~* 0. For large field intensities, the plateau 
extends linearly with the held amplitude in contrast with 
the case of high-order harmonic generation in ionizing 
systems, where the cut-off is proportional to the intensity. 

Figure [6] shows the results of Eqs. ( 22 ) and ( 24 1 for the 
case w/wo = 0.1 and different field amplitudes. Superim- 
posed to this, we plot the values of the harmonic number 
for the plateau's on-set and cut-off extracted from the 
numerical integration of Eqs. and For the latter 
case the plateau's on-set remains practically constant and 
equal to wo j however for larger photon frequencies, the 
departure between Eq. ( |22| ) and the low coupling limit wo 
becomes more apparent. For instance in the case plot- 
ted in figure [4] the plateau's onset frequency evaluated 



VI. PHYSICAL INTEPRETATION: THE 
ADIABATIC REGIME. 

The above formulas have a clear physical interpreta- 
tion if we consider the limit of small driving frequencies, 
w -C wo and w -C |%| [12j [13] . In this case we can con- 
sider the instantaneous eigenstates of the time dependent 
Hamiltonian as physically meaningful. The diagonaliza- 
tion of ([!]) and ^ gives the eigenenergies 



w ± 



Wn 



AF 2 {t) 



(25) 



with F(t) = Re {xexp(— iwi)}, oscillating harmonically. 
Figure [7] shows an schematic plot of these two eigenener- 
gies along a laser cycle. The instantaneous transition en- 
ergy is given by A + — A_ = ^/w 2 + 4F 2 (t) which oscillates 
between wo (at times when F(t) = 0) and \fwa + 4|x| 2 
(at the maximum amplitude of F(t)). These two val- 
ues coincide with the harmonic plateau on-set and cut- 
off at the limit w — > derived in the above section, see 
Eqs. (22 1 and (24). This suggests that the plateau struc- 



ture of the two-level dipole spectrum is originated by the 
transitions between the field dressed states. This mech- 
anism is different from the strong-field scenario, where 
the same structure is originated by the high energy ra- 
diation emitted from the recollision between an ionized 
electron with the parent atom. A fundamental difference 
between these two cases also appears in the identification 



8 



of the moment where the higher frequency harmonics are 
emitted. In the two-level case, Eq. (25), these harmon- 
ics correspond to the maximum of and, therefore, to 
the maxium of the driving field, while in real atoms the 
ionised electron recollides with the parent atom at times 
where the driving field is near zero. 



o< 











(a) ! (b) 









0.25 0.5 0.75 1 
Time (field period units) 



FIG. 7: Scheme of the instantaneous eigenenergies of the two 
level atom during one driving field period. In the adiabatic 
limit, the maximum (a) and minimum (b) level differences de- 
fine the energies of the plateau on-set and cut-off respectively. 



the electric dipole approximation for monochromatic 
fields. This new approach allows to perform a fundamen- 
tal decomposition of the complex dipole spectrum struc- 
ture into two families: the harmonic frequencies of the 
driving field and the field-induced nonlinear fluorescence 
frequencies. Moreover, it provides analytical expressions 
for the non-resonant Stark shift for the fluorescence fam- 
ily and of the on-set and cut-off limits for the harmonic 
family. Within this formalism, we predict the generation 
of pulses with lower duration than the period of the driv- 
ing field by selectively filtering out the frequencies of the 
emitted radiation, in close analogy with the attosecond 
pulse generation in the case of ionizing electromagnetic 
fields HI]. 

It is important to note that the presented formalism 
is an optimal tool to treat the burgeoning new physical 
scenarios where the two-level system approximation is 
applied, such as harmonic generation in molecules or 
dipolar two-level transitions in nuclei. Therefore, our 
approach is suitable to theoretically address near future 
experiments in the field. 
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